A Sequence of Quantum Gates 



(N 

o 

(N 



(N 



Yorick Hardy* and Willi-Hans Steeb^ 

* Department of Mathematical Sciences, 
University of South Africa, Pretoria, South Africa, 
e-mail: hardyy@unisa.ac.za 



f International School for Scientific Computing, 
^ ■ University of Johannesburg, Auckland Park 2006, South Africa, 

e-mail: steebwilli@gmail.com 



Abstract We study a sequence of quantum gates in finite-dimensional Hilbert 
spaces given by the normalized eigenvectors of the unitary operators. The corre- 
sponding sequence of the Hamilton operators is also given. From the Hamilton 
operators we construct another hierarchy of quantum gates via the Cayley trans- 



form. 
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J{~j: 1 Introduction 



In the finite dimensional Hilbert space C n quantum gates are given bynxn unitary 



matrices [H El El IH E] . For example for n = 2 we have the phase gate, NOT-gate 
and Hadamard gate 

Up= [l £)■ ^or=(j J), U " = ^(\ -1 

For the Mach-Zehnder interferometer these quantum gates play the central role 
[5]. The Pauli spin matrices a x , a y , o z with a x = Unot are also quantum gates. 
In higher dimensions we have the CNOT gate, Toffoli gate, Fredkin gate 

Ucnot = h ® cr x , U T = h@a x , U F = {l)@cr x @h 

where © denotes the direct sum. Teleportation in the Hilbert space C 8 is provided 
by the unitary matrix [5] 

U = U % UiUeU h UiUzU 2 U x 

where 



Ui = I 2 ®U H ® h, 



U2 — h <8> Ucnot, 



U 3 = Ucnot <8> h, U^ — Uh <8> h® h, 
U 5 = I 2 ® U C not, U 6 = I 2 ®h® U H , 
U 7 = h@U NO T@U NO T, U % = I 2 ® h®U H . 

The n x n unitary matrices form the compact Lie group U(n). If det(CZ) = 1 we 
have the Lie subgroup SU(n). The eigenvalues of a unitary matrix are of the form 
e i<t>j = i ; 2, . . . , n), where <pj G [0, 2tt), i.e. they are elements of the commutative 
Lie group £7(1). 

Beside the Hilbert space C" we also utilize the Hilbert space ofnxn matrices over 
C with the scalar product 

(A, B) :=tr(AB*). 

This implies a norm ||A|| 2 = ti(AA*) and a distance measure (Hilbert Schmidt 
distance) 

\\A - B\ = ti((A -B)(A- B)*) = tr(AA*) - tv(AB*) - tv(BA*) + tv(BB*) . 

In particular if we have two n x n unitary matrices, say U and V then 

\\U - V\\ =2n- tr(UV*) - tr(VU*) 

since tr(UU*) = n, tr(VV*) = n. Consider the two n x n unitary matrices U and 
e^U. Then 

\\U - e ^U\\ = 2n-2ncos0 

where <fi G [0, 27r). The norm takes a minimum for = 0. We consider two unitary 
matrices U and V equivalent if and only if V — e^U for <fi G [0, 2tt). This leads to 
equivalence classes. Thus for example a y and ia y belong to the same equivalence 
class. We define the distance measure [6] 

d(U,V):= min -^Llle^C/ - V II . 



We ask the following question. Given a unitary matrix Uq (not of the form e l ^I n , 
where I n is the identity matrix). We calculate the eigenvalues and correspond- 
ing normalized eigenvectors of U$. Using the normalized eigenvectors as columns 
we form a new unitary matrix U\. This process is repeated and thus provides 
a sequence of unitary matrices. Now every unitary matrix U can be written as 
U = exp(K), where K is a skew-hermitian matrix. K can be identified with a 
Hamilton operator (hermitian matrix H) via K = —iH. Thus we find the corre- 
sponding sequence of Hamilton operators H , Hi, .... 



From these hermitian matrices we can also form a new hierarchy of unitary matrices 
via the Cayley transform. The Cayley transform for annxn hermitian matrix H 
is given by 

V = (H-iI n )(H + iI n y l 

where V is a unitary matrix. Note that +1 cannot be an eigenvalue of the unitary 
matrix V. 



2 Sequence of Unitary Matrices 

Let { e 1>n , . . . , e n>n } denote the standard basis in C n . Let < denote an ordering 
on C n and let F n :U{n) —¥ U(n) be a function with the properties that F n (U)ej :n 
is an eigenvector of U for j = 1,2, ... ,n and the columns of F n (U) are in order 
with respect to <1, i.e. 

F n (U)e hn < F n (U)e 2 , n <■■•< F n {U)e^ n . 

Consider a sequence Uq (given), Ui, U2, ... of unitary matrices acting on C n defined 
by 

Uk+i ■= F n {U k ). 
We have the spectral decompositions 

U k := F n (U k ) diag (e^\e^\ e**.») F n (U k )* 

where e l6k < j is the j-th eigenvalue of U k with corresponding eigenvector F n {U k )ej^ n . 
It follows that the corresponding Hamilton operators are 

H k ■— F n {U k ) diag (—9k,u ~^k,2, ■ ■ ■ , —Ok,n) F n (U k )* 
and the Cayley transforms are given by 

V k := -F n (U k ) diag - — — S - — — S . . . , - — —2- F n (U k ) . 

\l — V kj l l — V k) 2 1 — v k) nJ 

Suppose that the sequence converges in the sense of the Hilbert Schmidt distance. 
Then 

lim U* k+1 U k = lim diag (e**-\ e ie ^, . . . , e"*-) F n (U k )* = I n 

k— >oo fe— ►oo 

from which follows 

lim U k = lim F n (U k ) = lim diag (e**- 1 , e^ 2 , . . . , e i6k ' n ) . 

k—>oo fc->oo k— yoo 



Thus for convergence we must have 




e 



for all «i, a 2 , ■ ■ ■ , a n G R. 

We note that even though the eigenvalues need never explicitly appear in the 
definition of F n , they still play a role in the convergence of the sequence. 

3 Examples 

To define the sequence we need 

1. the ordering <l, 

2. the function F n , 

3. and the initial unitary matrix U . 

We consider a few different U - For each of the sequences we will use the same 
ordering and the same definition of F n (dependent on n). 

Example <l. 

The ordering < we consider is 



x<y (K (x , y) x| > |< (Xiy) y|)v(|< (x>y) x| = |e; (Xjy) y| A arg(e; (Xjy) x) < arg« (Xjy) y)) 



w(x, y) := min {v G {1, . . . , n} : e*x ^ e*y} . 

Example F n . 

For each of the sequences we will use the same definition of F n (dependent on n). 
We choose F n (U) as follows. Let {x 1; x 2 , . . . , x„} be an orthonormal set of eigenvec- 
tors of U. The eigenvectors are determined as follows. Let I1 A denote the projection 
operator onto an eigenspace of U for the eigenvalue A. Then the Gram-Schmidt 
orthonormalization process applied to {n A e l n , n A e 2)n , . . .n A e n n }, in that order 
(and discarding zero vectors), yields an orthonormal basis {xj,x J+ i, . . .} for that 
eigenspace. Define 



where 



yfe := 



; l(x fc ),n X fc 
7(x fc ),n X fc 



Xfc, 



/(x) := min \y G {!,..., 



Now let a : {1, 2, . . . , n} — > {1, 2, . . . , n} be the permutation such that 



y<r(i) <y*(2) <---<ya(n)- 



Then we define 

n 

F n (U) := $>*Ci) e i,n- 

Example 1.1. 

As a first example we consider the U = U^ot = &x gate. This is a special case of 
example 1.2 below. The corresponding Hamilton operator H is 

Hq — - 



2 \-l 1 

and the unitary operator given by the Cayley transform of H is 



V n = 



1 in 



1 + iir \in 1 



with the eigenvalues +i and —%. The eigenvalues of £7 are +1 and —1 and the 
corresponding normalized eigenvectors are 



. -( 1 ■ 



This leads us to the Hadamard gate 

The eigenvalues of C/i are +1 and — 1 with the corresponding normalized eigenvec- 
tors 



J_ ( V^ + 2V2\ J_ ( ^4 - 2v/2 \ 
v/8 V \/4 - 2v/2 / ' ^1-^4 + 2^/ ' 



The corresponding Hamilton operator is 



4-2^ -v^ 



H =-- 

" 8 V -V8 4 + 2v/2 



with the Cayley transform of H 



V 1 



At the next step we find 



1 f-V2n-2i -y/2- 



7Y 



jj = J_ ( y/A + 2v/2 ^4-2^ 
2 V8 V \A - 2^ -^4 + 2^2 



and (a := 2 + V2, (3 := 2 - V2) 

rr Zll ( 1 (2 + v ^)/v^ \ 

2 8 + 4^^(2 + ^)/^ (4^ + ^2 + 6)//^ 
The distances between the unitary matrices are 



\\Uo-U4 = \/a-2V2, \\U 1 -U 2 \\ = \jA-\] A + 2^2-\] A- 2^2. 
From example 1.2 we find 

1 



lim U k = , 

k— >oo \ U — 1 



Example 1.2. 

Consider the unitary matrix 



a b 
b —a 



U = 

where a, b G R, b 7^ and a 2 + 6 2 = 1. If a > then we find 



sgn(6)6' -sgn(6)a' > /' V 2 ' V 2 ' 

and for a < we have 



f. = l *,„ u .J'.A * = /1 + a 



-sgn(6)6' sga(b)a' J ' ' V 2 ' V 2 ' 

while for a = we find that U\ is the Hadamard gate 

rr b> \ ,11 

Ul= {b> J' a: =7!' 6: =71- 

In each case, by construction, we have a' > and b' ^ 0. We also note that for 
a > we have a' > \a\ and b' < \b\. Since a', b' e (0, 1] we find that sequences 
beginning with 

a 6 



where a > 0, converge to a z . 



b —a 



Example 2. 

As our second example we consider the Pauli matrix 



U = cr 2 = 

with 



-i 

1 



Bo = -l(\ IV V 1 ' 1 » 



2 V ~i 1 / ' 1 + wr V _7r 1 

The eigenvalues of Z7o are 1 and —1 and we obtain the unitary matrix 

W i i 



The eigenvectors of XJ\ are 



y/2 \-l I 



V12 + 4V3 V -(! - 0(V3 + 1) J ' ^12-4^3 V (1 " 0(>/3 - 1) 
Thus 



Uo 



1 / 2^12-4^ 2 v / 12 _ f4V3 
4^1(1-0^24 + 873 (1-^^24-8^3 



Example 3. 

As last example we consider the 3x3 matrix 

l/y/2 l/y/2 



Un 



l/y/2 -l/v^ 



Then 



_ 'A-2\/2 — s/8 
H = ~\ 



8 



-y/E 4 + 2v/2 



x /-y/2n-2i -v^vr 

2(^ + ^ _y2vr V2tt - 2* 

The eigenvalues of ?7 are +1 and —1 with the corresponding normalized eigenvec- 
tors 



For the first step in the sequence we find 



1 / ^4 + 2^ \ 

U 1 = -=\ _0 v^8 • 

V8 \s/A-2y/2 -^4 + 2^ / 

4 Kronecker Product, Direct Sum and Star Prod- 
uct 

From the 2x2 case we can easily construct these sequences for the 4x4 case given 
by the direct sum, Kronecker product and star product. The star product of two 
2x2 matrices A and B is defined by 

/ an ai 2 \ 

A*B-= ° 611 612 ° 
6 2 i b 2 2 

\ a 2 i a 22 ) 

Note that if A and B are unitary, then A-kB is unitary We also have det(A*B) = 
det(A) det(B) and tr(A * B) = tr(A) + tr(B). 

In particular if F n distributes over any of these operations i.e. 

F mn {U A ®U B ) = F m (U A )®F n (U B ), 
F m+n (U A @U B ) = F m (U A )@F n (U B ), 
(U A *U B ) = F m (U A )*F n (U B ) 

for U A G U(m) and U B G U(n), then the sequence for Uq := C/0,0 <8> C/0,1 etc. 
is completely determined by the sequence for Uoo and Uoi- In other words if 
F m n{U A ®U B ) = F m (U A )®F n (U B ) then we construct the sequences for C/ ,o G C/(m) 
and C/0,1 G C/(n) 

^o,0) ^1,0 = ^1(^0,0)) ^2,0 = ^(^1,0), 
C/0,1, C/1,1 = F n (U ,i), C/2,1 = F n (C/i,i), . . . 

and 

C/ , C/i = F mn (U ) = C/1,0 <8) C/1,1, U 2 = F mn {U 2 ) = U 2fl ®U 2 ^ ... 



5 Conclusion 

We defined a sequence of unitary operators on the Hilbert space C n . The eigen- 
vectors of each unitary operator is used to construct the next. We found some 
necessary criteria for convergence. 

Some open questions remain. The role of the eigenvalues in the convergence of 
the sequence has not been determined. Due to the implicit role of the eigenvalues 
we have not examined the sequence of eigenvalues and their convergence properties. 

Obviously this construction can also be applied to other invertible matrices. For 
example the matrix (a G M) 



is an element of the non-compact Lie group 50(1,1) (Lorentz transformation). 
The matrix is also hermitian. Applying the construction described above we find 
the Hadamard gate Uh which is an element of the compact Lie group U(2) [7]. 
Since A (a) is hermitian we can apply the Cayley transform and find the quantum 

gate 



cosh(a) Vsinh(a) -i J 
with U(a — > oo) = o x = Unot- 
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